Saint- Venant's principle in dynamical porous 
thermoelastic media with memory for heat flux 

o 
o 

(N ■ Gerardo Iovane, Francesca Passarella 



(N 



o 



Department of Information Engineering and Applied Mathematics (DUMA), 
University of Salerno, 84084 Fisciano (Sa), Italy 



> 

ly-) ■ Abstract — In the present paper, we study a linear thermoelastic porous material 

with a constitutive equation for heat flux with memory. An approximated theory 

, of thermodynamics is presented for this model and a maximal pseudo free energy is 

determined. We use this energy to study the spatial behaviour of the thermodynamic 

processes in porous materials. We obtain the domain of influence theorem and establish 

the spatial decay estimates inside of the domain of influence. Further, we prove a 

uniqueness theorem valid for finite or infinite body. The body is free of any kind of a 

^ fi-i priori assumptions concerning the behaviour of solutions at infinity, 

-i— > 

1 Introduction 

^ ■ The Saint- Venant's principle has a central role in more theoretical and applied questions of 

elasticity. An important review of research on the spatial behaviour of solutions for statical 
and dynamical problems was given by Horgan and Knowles [1] and Horgan [2, 3]. Relevant 
information on the spatial behaviour of the solutions for the dynamical problems of elasticity 
are given by the domain of influence theorem as it is presented by Gurtin [4] . A further study 
in this connection was made by Chiri^a [5] for the linear theory of thermoelasticity, where a 
bounded solid is subjected to the action of nonzero boundary loads only on the end face. 

Recently, Fabrizio, Lazzari and Munoz Rivera [6] have studied a linear thermoelastic 
material which exhibits a constitutive equation for heat flux with memory. Some existence, 
uniqueness and asymptotic theorems have been established in connection with this approach. 
Further, Chirita and Lazzari [7] have studied the spatial behaviour of the thermoelastic 
processes. The hereditary effects are taken into account in a dissipative boundary condition 
by Ciarletta [8] and Bartilomo and Passarella [9]. 

*E-mail address: iovane@diima.unisa.it, passarella@diima.unisa.it 
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On the other hand, Iesan [10] has developed a linear theory of thermoelastic materials 
with voids, by generalizing some ideas of Cowin and Nunziato [11]. In this theory the 
bulk density is written as the product of two fields, the matrix material density field and 
the volume fraction field. This representation introduces an additional degree of kinematic 
freedom and it is compatible with the theory of granular materials developed by Goodman 
and Cowin [12]. The spatial and temporal behaviour in linear thermoelasticity of materials 
with voids has been studied by Chiri^a and Scalia in [13]. Further, it is shown by Iovane 
and Passsarella in [14] the spatial exponential decay for the boundary-final value problem 
associated to elastic porous materials with a memory effect for the intrinsic equilibrated 
body forces. 

In the present paper, we study a linear thermoelastic porous material which exhibits effects 
of fading memory in the constitutive equation for the heat flux. Such model exhibits finite 
speeds for the propagation of thermal disturbances and hence leads to accurate modeling of 
the transient thermal behaviour. By applying the thermodynamical laws, we can prove the 
existence of a pseudofree energy potential such that the internal dissipation vanishes despite 
of the dissipative character of the model. Then, we use this energy for studying the spatial 
behaviour of thermodynamical processes. 

The paper is organized as follows. In Section 2, we present the thermodynamical restric- 
tions imposed by the thermodynamical laws on the model in question; then, we introduce a 
maximal pseudofree energy. In Section 3, we deduce a result for describing the domain of 
influence by using an appropriate time-weighted surface measure. This domain at instant 
t G [0, T] is the set D^ t of all points of continuum for which r < (t, where r represents 
the distance from the bounded support Dt of the (external) given data in the time interval 
[0, T] and ( is a constant characteristic to the thermodynamic coefficients. We obtain, into 
the inner of the domain of influence, a spatial decay of exponential type. The decay rate, 
characterized by a factor independent of time, is exp(— — r), where a is a positive parameter 

independent of time and defining the time-weighted surface power measure. In Section 4, we 
prove that, into the inner of the domain of influence D^ t , an energetic measure associated 



2 Basic equations. Thermodynamic restrictions 

We consider a porous thermoelastic body that occupies a (regular) region B of the physical 
space IR 3 in an assigned reference configuration. By identifying IR 3 with the associated vector 
space, an orthonormal system of reference is introduced, so that vectors and tensors will have 
components denoted by Latin subscripts ranging over 1, 2, 3. The letters in boldface stand 
for the tensor L of an order p > 1, and L^-. . . k (p subscripts) are the components of the 
tensor L. Summation over repeated subscripts and other typical conventions for differential 
operations are implied, such as a superposed dot or a comma followed by a subscript to denote 
partial derivative with respect to time or with respect to the corresponding coordinate. 



with the thermodynamic process tends to zero with a decay rate equal to ( 1 — 4r 
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We denote U = {u, ip, 9}, where u is the displacement vector fields, ip is the change 
in volume fraction starting from the reference configuration, 9 is the temperature variation 
from the uniform reference temperature #o(> 0). Further, p is the bulk mass density and X 
is the equilibrated inertia in the reference state. 

We restrict our attention to the linear theory of thermoelasticity in which the effect of 
fading memory is contained in the thermal gradient V0. The local balance equations become 
[10] for the present problem 

S ji,j + Pfi = P^h 

hj,j + g + = pxfi, (i) 

—qjj + pr = pr, on B x (0, oo). 

In the previous equations: 

S and f are the stress tensor and body force, respectively; 

h, g and £ are the equilibrated stress vector, intrinsic and extrinsic equilibrated 
body force, respectively; 

r, q and r are the rate at which heat is absorbed for a unit of volume, the heat 
flux vector and the (extrinsic) heat supply, respectively. 

Let E be the strain field associated with u, y#' be the history up to time t for the thermal 
gradient and Xj9 be the integrated history of the thermal gradient, i. e. 

E ij = \( u ij + U 3,i)> ( 2 ) 



(3) 



O'jis) = 9j(t - s), 9%s) = f 9j(a)da, Vs e [0, +oo). 

Jt-s 

As it was shown by Ie§an in [10], the constitutive equations are 

Sij = Cij rs E rs + Dij S p yS + Bijip + Mij9, 

hi = D rsi E rs + A is p yS + bip + afl, 

(4) 

g = -B rs E rs - b s p tS -£<p- m0, 

pr = -9 M rs E rs - 9 a s p, s - 9 m<p + pc9, 

where the material coefficients of eqs. (4) are supposed continuous and bounded fields on B 
(the set B is the closure of B). Moreover, they satisfy the symmetry relations 

C - C 1 — C D — D A — A 

^ijrs ^rsij ^jirsi ^ijr ^jtri 13 311 

(5) 

Bij = Bji, Mij = Mji. 
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Further, we consider for the heat flux the following constitutive equation as suggested by 
Fabrizio, Lazzari and Munoz Rivera in [6] 

POO 

q t (t) = -0 O / K^e^ds, (6) 
Jo 

with the relaxation thermal conductivity K. We are assuming Kij are continuous and 
bounded fields on B; moreover, Kij(s) e L 2 (0, oo) and they satisfy the relations 

K ij (s) = K ji (s), K ij (oo) = 0. (7) 

By taking into account 0*-(O) = and 

^(s) = e d (t) ^(*) = *,-(«), (8) 

the constitutive equation (6) can be written as 

POD 

q t (t) = o / kM&jtfds. (9) 
Jo 

In what follows, we study the restriction imposed by the fundamental laws of thermodynam- 
ics in terms of cyclic processes. With this aim, we recall the laws of the thermodynamics as 
considered in [6, 15, 16]: 

First Law of Thermodynamics: For every cyclic process the following equality holds 
j jpr(t) + S rs (t)E rs (t) + h s (t)ip, a (t) - g(t)<p(t)} dt = 0. (10) 



Second Law of Thermodynamics: For every cyclic process the following equality holds 

1 j {pr(t)(e - e(t)) + q s (t)e,s(t)} dt < 0, (11) 

and the equality sign holds if and only if the process is reversible. 
The relations (6, 11) imply 

J (^J K^e^ds^j e ti (t)dt > o (12) 

for any cycle on [0, d), and the equality sign holds if and only if \j6 l is a constant history. 
If we put = h cos(u;t) + fcj sin(c^t) into the relation (12) with d = then we obtain 

POO POO 

' (kiKij(s)kj + kiKij(s)kj) cos(uis)ds + / (kiKij(s)kj — kiKij(s)kj) sm(u!s)ds > 0, (13) 

•/() 



where u is strictly positive and k iy hi are the continuous functions on B. 

Now, we introduce the Fourier transform, sine and cosine transforms of function / G 
L 2 (— oo, oo), i. e. 



(14) 



J — OO 

poo poo 

f(u) = / f(0 sinK)^, f(u) = / /(0 cosK)^- 
Jo JO 

We remark that the Fourier inversion formula imply 

1(0 = ^11°° f(OM^)doo; (15) 
and the Plancherel's theorem of the Fourier transform gives 

/OO / 1 poo 

x mwm = v 2^7 oo f*® 9 **®**' (16) 

with f,g<E L 2 (-oo, oo) and with g* F is the complex conjugate of gp. 
We can easily prove by eqs. (8, 14) that 

9f = -u9f, 9f = ^(t)+^f. (17) 

By putting \ = ki=ki in the relation (13) as Chirita and Lazzari in [7], we obtain 

kiKfjkj > 0, Vk = (fci, k 2 , k 3 ) ^ 0, Vu > 0. (18) 

Thus, the relations (17, 18) imply that K s is a negative definite tensor. 
By using eqs. (9, 15, 16), we have 

= \J\e Q J™ k^uWfi^du;, (19) 



and 

K l3 (0 = \l-l mu)an^u)du. (20) 



If we integrate eq. (20) with the respect £ and we take into account the Riemann Lebesgue 
lemma and the hypotheses K^oo) = 0, then it follows 

/2~ f°° 1 
- y o -Kf^du. (2i) 
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The relations (7, 21) and the negative definiteness of the tensor K 5 imply that K(0) is a 
symmetric and positive definite tensor, i.e. 

kiKij^kj > 0, Vk = (h,k 2 , h) ^ 0, Vcj > 0. (22) 

Consequently, it follows that 

Kij(0)kikj < K M kih, Vk ={fci, k 2 , k 3 }, (23) 

where Km(x) > is the largest characteristic eigenvalue of K(x, 0). 

In what follows, we will assume that the bulk mass density p, the equilibrated inertia x 
and the constant heat c are strictly positive, continuous and bounded fields on B, so that 

< p = inf p, < xo = inf x, < c = inf c. (24) 
xeB xeB xeB 

Moreover, we suppose that K(x, 0) is continuous and bounded on B. By setting 

_ 9qK m 

K = sup , (25) 

xeB c 

we obtain 

i^(x, 0)*^ < Ko-^-hh Vk = {fci, fc 2j k 3 } ± 0. (26) 

Now, we introduce the vector space V of all vector fields of the form 

F = {Fij, y/x^i, ip}, with Fij = Fji. 

For any F, F G V, we define the inner product and the magnitude by 
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F tj + x*m + i>\ (27) 



F F = F^ + x^i + W>, 
where F = {F^, y/x^i,i>}- For given F G V, we introduce the vector field S(F) G V as 

S(F) = |^,(F), (^( F )) > -5(F)} , (28) 

where 

hi(F) = D rsi F rs + AijTTj + bit/;, (29) 

g(F) = -ByFy - bi-Ki - 
and the coefficients are defined in eqs. (4). They obey to the symmetry relations (5). 
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As in [13, 14], we consider the following bilinear form 

+ £V0 + AjKiiTj + B ij (F ij if) + FyV>) + 
+D ijs (F ij 7r s + FijTCs) + biiipWi + ipiXi). 

From (29, 30) we have 

2^(F, F) = ~S l3 {F)F l3 + hi(F)7ti - g(F)$. 



(30) 



(31) 



Throughout this paper, we will assume that the quadratic form W associated to T is positive 
definite, so that 

2W(F) = ZF(F, F) < n M (F ij F..+ X mm + ^ 2 ) , 



where /im > 0. In particular, in the case E = {Eij, y/x¥,i-, V 5 }? from (4, 29) we get 

Sij = S^E) + MijO, h i = h i (F) + a i e, g = g(E) — m9. 
Taking into account eqs. (29, 30, 32, 33), we deduce 

W* = W(E) = i[(% - Mij6)Eij + (h, t - OiO)^ -(g + m0)<p], 



(32) 
(33) 



W = SijEij + hiip ti - gip + — —6, 

"o 



and 



2W* < ii m {E 1] E ij + X <p,i<p,i + p 2 ) . 
By setting F = S(E) in the relations (31,32), we have 



2W(S{F)) < fi M 



4(E)4(E) + -h^kiF) + g 2 (E) 

X 



(34) 
(35) 

(36) 



Therefore, we conclude thanks to (31, 36) and thanks to Cauchy-Schwarz's inequality that 



^(E)^(E) + -^(E)^(E) + g 2 (E) 



X 



2F(E, S(E)) < AW*W(S(E)) < 



< 2W*hm 



%(E)%(E) + -^(E)^(E) + g 2 (E) 

X 



Consequently, it follows 



_ _ 2 



S(E) =,%(E),%(E) + -fn(B)hi(E) + g 2 (E) < 2pt M W*. 

X 

Thanks to the inequality for second-order tensors L and G 

(Lij + Gij)(Lij + dj) < (1 + e)LijLij + (1 + ^G^, Ve > 0, 



(37) 



and thanks to eqs. (33, 37), we get 

SijSij H hihi + 
X 









<(! + *) 


S(E) 





X 



< (l + e )2^* + (l + i)M^^ 2 , Ve>0, 



(38) 



with 

/i = sup (Km, M = sup — I MijMij H — a^di + m 2 J . (39) 

The Laws of thermodynamics imply the existence of the thermodynamical internal potential 
energy e and the entropy r] such that 

pe(t) = pr(t) + %(t)^(t) + / ij (t)^(t)-( ? (t)v9(t), (40) 

PV(t) > ^{pr(t)(9 -e(t)+ qj (t)9,(t)}. (41) 

We introduce the pseudofree energy potential 

p*(t)=pe{t)-p9 ri(t). (42) 
Moreover, we define the maximal pseudofree potential energy ty M such that 

P i> M {t) = + s^mm + hMMt) - 9 {mt) - q jMM. m 

Then, the relations (40-43) imply 

p*(t) < p* M (t). (44) 
By taking into account the relations (4, 19, 34 2 ), a maximal pseudofree energy potential is 

p* M (t) = w\t) + p -^- -lJ\{k^,u)ef^,u J )9 t J^,uj)+ (45) 
+ii:g(x, w )flf(x, w )^(x, w )}d w . 

Since W is a positive definite quadratic form, p, \ an d c are strictly positive, K 5 is the 
negative definite tensor, we deduce that the functional ^ m is a norm. As in [7], we can 
prove by eq. (19) that 

POO 

|q(t)| 2 = ft (t)/r*(x, W )^(x,u;)dw< (46) 

^ Jo 

(2 /"°° 1 \ !/ 2 / o roo _ _ \ 1/2 

7 -«(*)*(*)^HdwJ (-- jf ( W )^f(a;)du;J . 

The relations (21, 25, 26, 45, 46) imply 

|q(t)| 2 < K 9 c ^ M (t) - 2W*(t) - . (47) 
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3 A time-weighted surface power measure 



Throughout this work by an admissible process we mean an ordered array 7r = [u, E, S, 95, 7, 

h, 9, k, q] with the following properties 

i. u u <p G C 2 ' 2 (B x [0, +00)), 9 G C^{B x [0, +00)); 

ii. En = Eji, ^=ip :i , k t = 9\ { G C°'\B x [0, +00)); 

iii. Sij = Sji, K, qi G C 1 >°(B x [0, +00)), 

and which meets the equations of motion (1), the geometrical equations (2, 3), the consti- 
tutive equations (4, 6) and the following initial conditions 

u t (0) = <, «i(0)=«?, y?(0) = y?°, 0(0) = 0°, 0(0) =tf . (48) 

Let s, ft, and q be the surface tractions, the surface equilibrated stress and the heat flux, so 

Si(t) = Sjii^rij, h(t) = hj(t)rij, q(t) = qj{t)rij. (49) 

We denote by T = [/,, Sj, /, ft, <p°, 0°, r, g, O) 0°-] the external data and assume 

that all functions are prescribed continuous functions. 

Now, we introduce the support D T of the external data T and the body supplies on the 
time interval [0, T], i. e. the set of all x G B such that: 

(1) if x G B, then 

^ or u? ^ or V or V or O 7^ 

or 

0j( s ) ^ for some s G (-00, 0] (50) 

or 

fi(s)^0 or £(s)^0 or r(s) ^ for some s G [0 , T] ; 

(2) if x G dB, then 

Si(s)Mi(s) 7^ or h(s)ip(s) 7^ or q(s)9(s) 7^ for some s G [0,T]. 

In what follows, we will assume that Dt is a bounded set. 

We introduce a nonempty set D T such that Z) T C D T C -B and 

(a) if D T nB 7^ 0, then we choose D T to be the smallest bounded regular region in B that 
includes Dt; in particular, we set D T = Dt if Dt it is a regular region; 

(b) if 7^ Dt C <9£>, then we choose D T to be the smallest regular subsurface of dB that 
includes Dt'-, in particular, we set D T = Dt if Dt is a regular subsurface of dB; 

(c) if D T = 0, then we choose D T to be an arbitrary nonempty regular subsurface of dB. 
On this basis, we introduce the set D r , by 

D r = {xGB: D T fl E(r) 7^ 0}, 



9 



where E(r) is the open ball with radius r and center at x. 

Further, we shall use the notation B r for the part of B contained in B \ D r and we set 
B(ri, r 2 ) = B r , 2 \ B ri , r\ > r2- We denote by S r the subsurface of dB r contained into the 
inner of B and whose outward unit normal vector is forwarded to the exterior of D r . We 
can observe that the data are null on B r , S r . 

We define the following time-weighted surface power function J(r, t) 

I(r,t) = - f [ e-^is^u^s) + h(s)<p(s) - ^-q(s)9(s)]dads. (51) 

J0 J S r ^0 

for a fixed positive parameter a and for any r > 0, t e [0, T] 

The following theorems establish a set of properties for the surface power function /. These 
theorems will be useful in the study of the spatial behaviour of the thermoelastic processes. 

Lemma 1: Let n be a thermoelastic process and Dt be the bounded support of the external 
data T on the time interval [0, T}. Moreover, let I(r, t) be the time-weighted surface power 
function associated with tt and JC be the kinetic energy defined by 

K=^{pu lUl +pxf). (52) 

Under the hypotheses of the Section 2, it follows 
(I) for <r 2 < ri 



I(r u t) - I(r 2 ,t) = - f e-^[IC(t)+p^ M (t)}dv- 

t JB(ri,rj) 

-a I I e-° s [/C(s) + P^m(s)} dv ds; 

JO JB(r-i.ro) 



(53) 



'0 JB(n,r 2 ) 

(II) I(r, t) is a continuous differentiate function and 



^Hr,t) = -J^ e- CTt [/C(t)+p* M (t)] da- (54) 

-a [ [ e~ us [K{s) + p^m(s)] dads, 
Jo J s r 



and 

d 



I(r,t) = - f e-^ls^u^ + hit^-^q^e^da-, (55) 

J Sf 

(III) for each fixed t e [0, T], I(r,t) is a non-increasing function with respect to r. 



10 



Proof. By using the divergence theorem, by eqs. (1-4, 5, 49, 45, 29) and thanks to the 
definitions of D T , B r , S r , we have for < r 2 < r\ 

I(n,t) - I(r 2 ,t) = - [ [ e- as [s l {s)u l {s) + h{sMs)-^q{s)e{s)]dads = 

JO JdB(n,r 2 ) y 



* f 1 

/ e-^Sjij^iiiis) + h jJ (s)<p(s)- T 9(8)q iti (s) + S ij ( S )E ij (s) + 

JB( ri ,r 2 ) Un 



+hj(s)(pj(s) - —q j {s)e d (s)}dvds 
Vo 



// 

Jo Je 



d 



B(r 1: r 2 ) 



e-™— [)C(s) + pV M (s)] dvds. 



(56) 



From (56) we get eq. (53). 

The relation (53) implies eq. (54), taking into account the concept of derivative as incre- 
mental ratio and performing the limit r 2 — > ri; while, the eq. (55) comes from the definition 
of/(r, t). 

The property (III) is obtained from (I) and p, \ >0, and by considering that \lV is a 
norm. ■ 

Lemma 2: Under hypotheses of Lemma 1, the surface power function I(r, t) associated 
with n satisfies the following first-order differential inequalities 



d_ 

dt 



I(r,t) 



+ C|;/M) < 0, 



(57) 



^\I(r,t) \ +frI{r, t) < 0, 



where 



and 



( 1+foV 
Po 



1 (Kp + M) 
1 + £o= 2 + ~^n— + 



1 (K + M) y M 

2 '" 2p J p' 



(58) 



(59) 



(60) 



Proof. By using the Schwarz's inequality and the arithmetic-geometric mean inequality 
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and the relations (24, 38) and (47), it follows that 



Si (t)ui(t) + h(t)(p(t) - ^-q(t)6(t) 



< 



—pUi(t)ui(t) + — Sij(t)Sij(t)+ 
Po £i 



Po £i X Po y o £2 ctf 

< ^-/C(t) + \ \y ((1 + e)2^*(t) + (1 + \)M^$- ] + 

PO "0 £2 V PO 

< £L/C(t) + p M/ M (t) + f— (1 + e)p - - K ) W*{t) + 

Po L £2 \£i £2 / 

2 \ £ i e Po £2 / Oo 



< 



(61) 



< 



for every si > 0, £2 > 0, e > 0. We choose £, £1, £2 such that 

£l ^0 1 ^ . v ^0 n 1 

— = , — (l + e)p = 0, - 

Po £2 £1 £2 £1 



-(1 



-)M + — 
e Po 



^0 

£2 



0; 



consequently, we obtain 



£1 = PoC, 



£2 



Ko 

c ' 



where £ is given by (59) and £ is given by (60). In other words, £ is the positive root of 
the algebraic equation 



K + M M 

£ -1 + 

H p 



0. 



(62) 



Thus, the relations (51, 55, 61) imply 



I(r,t) 



< (f e- at [IC(t) + p^ M (t)}da, 

JSr 



a\I(r,t)\ < a ( / e~ as [K{s] + pV M (s)] dads. 

JO JSr 



(63) 
(64) 



By eqs. (54, 63, 64), the inequalities (57, 58) are reached. ■ 

Lemma 3: Under hypotheses of Lemma 1, the surface power function I(r, t) associated 
with 7r is equal to 

I(r, t)= f e-^[IC(t)+p^ M (t)] dv+ 



+a [ [ e' as [/C(s) + P^m(s)} dvds > 0, 

J0 J B r 



(65) 
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for r > 0,t G [0, T]. 

Proof. If B is a bounded body, then the variable r ranges on [0, L], where 

L = max xe jj{mffl y ^, |x — y|} < oo. 

From the definition of D T and 7(r, t), we have 7(7, t) = 0, so that we obtain eq. (65) 
with the help of eq. (53). 

On the other hand, if B is an unbounded body, then, the variable r ranges on [0, oo). 
Fixed a pair (r , to) in the plane (r, t), such that r > (\t — t \, t G [0, T], the inequality 
(57) implies that 

£[I(r + C(t-t ),t)] < 0; (66) 

l[I(ro-((t-t ),t)} > 0. (67) 

Therefore, we have 

I(r + t(t-t ),t) < I(r + ((t -t ),t) with t>t\ (68) 

I{r -((t -t ),t) > I(r -({t-t ),t). witht>t. (69) 

For t = t and t' = 0, the relations (51, 68, 69) imply 

I(r ,to) < 7(r -C*o,0) = 0, (70) 

= 7(r + C*o,0) >J(r ,t ), (71) 

and so 

7(r ,t ) = 0. (72) 
For r — > oo and, consequently, for any t G [0, T], eq. (72) becomes 

7(oo,t ) = lim 7(r ,t ) = 0. 

ro^oo 

From eq. (53) and for r\ — > oo we deduce 

I(r,t) = I(r,t) - 7(oo,t) = ^ e - ffl [IC(t)+^ M (t)] 



(73) 



By the properties of the surface power function 7(r, t) established in Lemma 1-3, we obtain 
a complete description of the spatial behaviour of the elastic process in question outside of 
the support of the data. 
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Theorem 1: Under hypotheses of Lemma 1, for each fixed t G [0, T] we have the following 
results: 

1) Domain of influence 

I(r, t)= 0, forr>(t. (74) 

2) Spatial decay 

I(r, t) < e-f r I(0, t), forr<(t. (75) 

Proof. By putting r = 0, t = in (66), it follows 

dJ(CM) 



dt 

and so 



<0, 



7(CM)</(0,0) = 0. (76) 
From Lemma 1 and Lemma 3, we have 

0<I(r,t)<I((t,t), Vr>Ct. (77) 

Therefore, by the inequalities (65, 76, 77) we obtain the equation (74). 
On the other hand, the inequalities (58, 65) imply 

A[ e ! r J(r,*)] <Q, forr<Ct, (78) 



so that we arrive to (75). ■ 

We remark that the relations (45, 65, 74) give 



= I(r,t)= I e~ at {)C(t) + [W*(t) + 



pc6 2 {t) 
29 



1 f°° 

; " a;{XgH^(a;)^(a;) + Xg(a;)^H^(a;)}da;]}^, for r > (t. 



7T 



Taking into account that p, x an d c are strictly positive functions, that If is a positive 
definite quadratic form and that K 5 is negative definite tensor, we have 

Ui = 0, = 0, 6 = on B r x [0, T]. (79) 

Since the external data are null on B r , the relation (79) yields 

Ui = 0, y? = 0, 6 = on B\D r x [0, T] for r > (t. (80) 
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Now, we establish the domain of influence of the external given data at time T, according 
to Gurtin [4]. In fact, we show, by putting t — T and r = (T in relation (80), that on [0,T] 
the external given data have no effect on points outside of D^t- 

Theorem 2: Under hypotheses of Lemma 1, it follows 

Ui = 0, y? = 0, # = onB\%x[0,T]. (81) 

As an immediate consequence of Theorem 1 , we show the following uniqueness result valid 
for a bounded or unbounded body: 

Theorem 3: Under hypotheses of Lemma 1, there exists at most one solution for the 
boundary-initial-value problem. 

Proof. Thanks to the linearity of the problem, we have only to show that the null 
data imply null solution. Let U = jttj,</5, #j be a solution of the problem (1, 6, 49, 48 ^ 

corresponding to null data. In this case, for each T G (0, +oo) we have D T = and 
I(r, t) = 0. Then, we can conclude that 

u, L = 0, <^ = 0, 6 = on B x [0,+oo). 



4 An alternative surface measure 



In the above section, we have assumed that the parameter a is strictly positive. We can 
avoid this restriction by considering a measure of the type 



P(r. t) = - I [ [Si(s)ui(s) + h(s)^{s) - ^q{s)6(s)]dads, 



JS. 



0, 







(82) 



for r > 0, t G [0, T]. Following the procedure developed in the previous sections, we can 
easily prove that P(r, t) satisfies 



§~ r P (r,t) < 0, 



(83) 



Moreover, we can deduce that P(r, t) is equal to the total energy associated to 7r on B T) i.e. 



P(r,t)= I [Kl(t) + pV M (t)] dv>0, 

I B r 



and 



P(r,t) = 0, for r > (t. 

On the basis of these results, we establish the following spatial decay result 



(84) 
(85) 
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Theorem 4: Let n be a thermoelastic process, D T be the bounded support of the external 
data T on the time interval [0, T] and P(r, t) be the surface measure associated with n. 
Under hypotheses of Section 2, for each fixed t G [0, T] we have 

Q(r,t) < (l-^Q(0,t), forr<(t, (86) 

where t 

Q(r,t)= P(r,a)da = / [/C(a) + p# M (a)] dvda. (87) 

JO J0 J B r 

Proof. By (85, 87) we have 

ft' 



Q(r',t') = J rf P{r',a)da, 



for any choose r' > 0, t' e (0, T] so that r' < (t'. Now, we consider the following changed 
variable 



r'\ r' 



clearly, a > s. We get 



Q(r', = (l - ^) jf* P(r', (i " 7^) s + (89) 



On the other hand, as in the Lemma 3, we can prove that 

d 



dt P(r + ((t-t ),t) <0. (90) 



(91) 



Thus, by putting r = r', t — a into (90), we reach the result 

^,<,)<P(/ + «.-„),.) = p(^,.)^h„>.. 

The first inequality of (83) implies that P(r, t) is non-increasing function of r, such that 

P (^-,s^ < P(0,s) with s G [0,t'\. (92) 
For r' < (t f and t' G [0, T], we deduce by the inequality (91, 92) that 

p ( r/ '( i -^) s+ 9 - p(M ' } - (93) 

Finally, the relations (87, 89, 93) yield to the inequality (86). ■ 

However, we have to outline here that the decay rate characteristic to the estimate (86) 
is lower than the established one in the above section. 
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